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Local spin resonance and spin-Peierls-like phase transition 
in a geometrically frustrated antiferromagnet 
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Inelastic magnetic neutron scattering reveals a localized spin resonance at 4.5 meV in the ordered 
phase of the geometrically frustrated cubic antiferromagnet ZnCr2 04. The resonance develops 
abruptly from quantum critical fluctuations upon cooling through a first order transition to a co- 
planar antiferromagnet at T c = 12.5(5) K. We argue that this transition is a three dimensional 
analogue of the spin-Peierls transition. 



PACS numbers: 76.50.+g, 75.40. Gb, 75.50.Ee 

It appears that antiferromagnetically interacting 
Heisenberg spins on the vertices of a lattice of corner- 
sharing tetrahedra cannot order [Q-Q. Characterized by 
weak connectivity and frustrated interactions || , this so- 
called pyrochlore antiferromagnet has no phase transition 
for spin S — oo |^|,|3| and forms a quantum spin liquid 
with low lying singlet excitations for S = 1/2 jj]. Re- 
cent experimental |j]|7| and theoretical |^,^| work how- 
ever indicates that small perturbations away from the 
ideal model can induce low temperature phase transi- 
tions. In the quantum critical spin- 1/2 antifcrromagnetic 
chain sensitivity to small symmetry breaking perturba- 
tions leads to the spin-Peierls transition: A lattice distor- 
tion that lowers the energy of the spin system by inducing 
a gap in the magnetic excitation spectrum. In this letter 
we report a similar phenomenon for spins on a lattice of 
corner-sharing tetrahedra. 

We examined the spinel antiferromagnet (AFM) 
ZnCr 2 04 using neutron scattering. Cubic and paramag- 
netic at high temperatures, the material undergoes a first 
order phase transition at T = 12.5 K into a tetragonal 
phase with Neel order. Our measurements and analysis 
indicate that lattice energy is expended at the phase tran- 
sition to break the symmetry of frustrated magnetic in- 
teractions and thereby select an ordered magnetic phase 
from a manifold of degenerate states. A local spin reso- 
nance develops abruptly below the transition. This res- 
onance is a manifestation of Q-space degeneracy in the 
ordered phase of a frustrated magnet. Just as the de- 
velopment of a singlet-triplet gap drives the spin Peierls 
transition, the fact that ordering in a frustrated magnet 
can push low energy spectral weight into a finite energy 
resonance plays an important role at the transition in 
ZnCr204. 

Cr 3+ is the source of magnetism in ZnCr204. Approx- 
imately octahedrally coordinated, the unfilled 3d 3 shells 
of these atoms form isotropic 5 = 3/2 degrees of freedom 
on a lattice of corner-sharing tetrahedra. Assuming that 



only nearest neighbor interactions are appreciable yields 
an estimate of J = (31cb@cw / zS(S + 1)) = —4.5 meV 
for the exchange constant. Here Qcw ~ —390 K is the 
Curie- Weiss temperature derived from high temperature 
susceptibility data ||] and z — 6 is the nearest neigh- 
bor coordination number. Experiments on Cr 3+ pairs in 
ZnGa2 04 0,|ll| give values for J ranging from —2.8 meV 
to —4.0 meV and also provide evidence for biquadratic 
exchange H 2 = j(Si ■ S 2 ) 2 with j = -0.21(4) meV. 

A 25 g powder sample was prepared by solid state 
reaction between stoichiometric amounts of Cr2 03 and 
ZnO in air. Rietveld analysis of neutron powder diffrac- 
tion data from the NIST BT1 diffractometer shows that 
ZnCr204 in the spinel structure (space-group Fd3m, 
a = 8.31273A for T = 15 K) is the majority phase with 
a minority phase of 1% f. u. un-reacted Cr203. Elas- 
tic and inelastic neutron scattering measurements were 
performed at NIST on the cold neutron triple-axis spec- 
trometer SPINS. A vertically focusing Pyrolytic Graphite 
(002) monochromator (PG(002)) extracted a monochro- 
matic beam with energy 2.5 meV< Ei < 14 meV from 
a 58 Ni coated cold neutron guide. The detection sys- 
tem consisted of a 20 cm long polycrystalline BeO fil- 
ter cooled to 77 K followed by a 23 cmx 15 cm flat 
PG(002) analyzer 92 cm from the sample, then a 80' 
radial collimator, and an area sensitive detector. The en- 
ergy range detected was 2.6 meV< Ef < 3.7 meV with 
Full Width at Half Maximum (FWHM) energy resolu- 
tion 0.1 meV< AE < 0.15 meV and angular resolution 
A29 w 50'. The absolute efficiency of the instrument 
was measured using incoherent elastic scattering from 
vanadium and nuclear Bragg peaks from ZnCr204. The 
corresponding normalization factor was applied to back- 
ground subtracted data to obtain measurements of the 
normalized magnetic neutron scattering intensity |L3] 
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Here F(Q) is the magnetic form factor for Cr 3+ O] and 



K = 0.48(5) A" 1 = 0.64(6)a* indicating that the excited 
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Fig. 1. Color maps of the magnetic neutron scattering inten- 
sity versus wave vector and energy transfer at three temper- 
atures spanning the phase transition at T c = 12.5(5) K. 

5 Q/3 (Q,w) is the scattering function fll2| . 

Fig. 1 provides an overview of our data in the form 
of color images of I(Q,u>) at three temperatures. In the 
paramagnetic phase, Figs. 1 (a) and (b) provide evidence 
for quantum critical fluctuations of small AFM clusters 
most likely antiferromagnetically correlated tetrahedra. 
The data closely resemble those obtained in similar ex- 
periments on other frustrated AFM's @@. For T < T c 
however, the low energy spectral weight concentrates into 
a sharp constant-energy mode centered at %uj = 4.5 
meV« | J | >> ksT c . Above the resonance there is an 
additional band of intensity centered at fuv w 9 meV 
(Fig. 2 (a)). Careful inspection of data below the res- 
onance also reveals weak dispersing streaks emanating 
from AFM Bragg points. 

Focusing first on the ordered phase, Fig. 2 (a) shows 
the Q— integrated magnetic scattering cross section at 
T = 1.7 K. With a resolution corrected FWHM of only 
0.8 meV the Tuo = 4.5 meV peak is evidence for a 
near dispersionless excitation. The spectral weight is 

*<f!Z V evT, aa S aa (")d" = 0-59(1) /Cr which is 22% 
of the total fluctuating moment ((5(5 + 1) - |(S)| 2 ) = 
2.65(5)/CV) and corresponds to 16% of the total mag- 
netic scattering cross section. Fig. 2 (b) shows the Q- 
dependence of the energy integrated intensity of this 
mode. The data exhibit a broad peak centered at 
Qo = 1-5 A" 1 with a Half Width at Half Maximum 
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Fig. 2. Integrated magnetic scattering intensities at T = 
1.5 K derived from Fig. 1 (c). (a) oj-dependence of the Q- 
integrated intensity: I(u) = J Q 2 dQI(Q,cj)/ J Q 2 dQ. The 
horizontal bar shows the instrumental energy resolution, (b) 
Closed symbols show the Q-dependence of the ^-integrated 
resonance intensity: I(Q) = h J^'° I(Q,uj)duj. Open cir- 
cles show data integrated from 1 meV to 3 meV. The solid line 
shows the elastic scattering cross section scaled by a factor 10. 
Shaded peaks are magnetic Bragg peaks. 

state involves AFM correlated nearest neighbor spins. 
For comparison the dashed line shows the powder- 
averaged magnetic neutron scattering intensity for an 
isolated spin dimer at the nearest neighbor separation 
7"o = 2.939 A jl5|]. The spin pair model produces a 
broader peak than the experiment indicating that the res- 
onating spin cluster in ZnCr204 is more complex. Fig. 2 
(b) also shows the wave vector dependence of inelastic 
scattering below the resonance integrated over energy 
from 1 meV to 3 meV (open symbols). Weak, non- 
resolution-limited peaks are visible and their locations 
coincide with AFM Bragg peaks (shaded). Excitations 
below the resonance are also apparent in the wave vector 
integrated data of Fig. 2 (a) where the intensity increases 
in proportion to energy for Tim < 3.5 meV (the low en- 
ergy upturn is incompletely resolved elastic scattering). 
Both features are consistent with neutron scattering from 
spin waves in a three dimensional AFM powder with a 
spin gap A < 1.5 meV. In particular the ratio of elastic 
to inelastic scattering is consistent with estimates based 
on spin wave theory in the long wavelength limit. From 
the width of the peaks we estimate a spin wave veloc- 
ity v — 18(2) meVA. This number is much less than 
the spin wave velocity for a bi-partite simple cubic AFM 
with J = -2.8 meV, v = (2/V3)z\J\Sa = 239 mcVA 
but only slightly larger than the spin wave velocity for 
a cubic AFM with the critical temperature of ZnC^C^: 
V = 2V3k B T c a/(S + 1) = 12.4 meVA. 

Spin wave theory provides a useful starting point for 
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understanding the resonance. Geometrical 
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Fig. 3. (a)-(c) x(Qo-! u ) at Qo = 1.5A -1 derived from mag- 
netic neutron scattering data via the fluctuation dissipation 
theorem. The solid lines are fits as described in the text, (d)- 
(e) Temperature dependence of the relaxation rate, Tq (T), 
and the inverse susceptibility, Xq^ (T) derived from the fits. 

frustration leads to constant energy surfaces or volumes 
for spin wave dispersion relations in reciprocal space. 
Such Q-space "degeneracy" in turn yields pronounced 
van-Hove singularities in wave vector averaged spectra. 
Reimers et al |l6) showed that the pyrochlore AFM has 
two degenerate modes for any Q in the Brillouin zone. 
Wave vector independent excitations also exist for the 
kagome AFM jllj and these have a real space interpreta- 
tion in terms of the so-called weather-vane modes |Q . A 
real-space interpretation has yet to be found for disper- 
sionless excitations in the pyrochlore lattice. The broad 
peak in Fig. 2(b) indicates that they are highly localized 
in the ordered phase of ZnC^O^. 

Turning now to excitations in the paramagnetic phase, 
Fig. 3 (a)-(c) show the imaginary part of the spin sus- 
ceptibility, x"{Qj w ) f° r several temperatures larger than 
T c - x'^Qi^) was derived from inelastic neutron scatter- 
ing data at Qq = 1.5 A -1 via the fluctuation dissipation 
theorem: x"{Q^) = (ff^s) 2 ^! - exp(-/?w))S(Q, w). 
From the spectra we derived a temperature dependent 
spin relaxation rate, Tq, and a static staggered suscep- 
tibility, xq, by fitting to the following phenomenological 
response function: x"(Q> u ) — XqFqu/(uj 2 + Tq). 

Figs. 3 (d) and (e) show the corresponding tem- 
perature dependent parameters. A power-law in T 
describes the temperature dependence of Fq (T) — 
C ■ kBTiT/e) 01 - 1 for T < 150 K while XQ ( T ) 
can be described by a Curie- Weiss law: XQo(T) = 
(/Xq o /3£;b#)1/(1 + T/9) in the entire temperature range. 
The best fit solid lines correspond to a = 0.81(4), 
C = 0.6(1), n Qo = 4.0(l)/i B , and 9 = 8.8(4) K. Though 
we have plotted and analyzed data for Qo = 1-5 A -1 we 
found similar results for T(T) in analysis of wave- vector 



integrated data which probe the local spin susceptibility. 
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Fig. 4. (a) Color image of inelastic neutron scatter- 
ing for Q — 1.5 A -1 . (b) T— dependence of magnetic 
Bragg scattering from a powder (blue squares), <r m = 
pffia J I(Q,uj)4:TYQ 2 dQdhuj where v m is the volume per Cr i+ 
ion, and of lattice strain along a and c (red circles) measured 
by single crystal neutron diffraction. 

At a quantum critical point, kgT is the only low en- 
ergy scale for local response functions. If the lorentzian 
form for ^"(Q,o;) describes the spectrum, then a must 
equal one at the critical point. This was the exponent 
found in Monte Carlo simulations of classical spins on a 
pyrochlore lattice ||. The deviation of a from unity in 
ZnCr204 is perhaps not surprising given that the mate- 
rial does exhibit a magnetic phase transition. However, 
the fact that Tq tends to zero as T — > rather than at 
T c indicates that the system may actually be approach- 
ing a quantum disordered phase with a gap ksO = 0.75 
meV before being interrupted by a first order transition 
to an unrelated competing phase. This idea is consis- 
tent with Fig. 4 which compares the temperature depen- 
dent lattice strain and magnetic Bragg peak intensity 
(frame (b)) to the inelastic neutron scattering spectrum 
at Qo = 1-5 A -1 (frame (a)). Magnetic Bragg peaks, a 
tetragonal lattice distortion, and the spin resonance all 
appear abruptly, and without conventional critical fluc- 
tuations at T c . 

Theoretical work has shown that magnetic order can 
not develop in an isotropic spin pyrochlore AFM ||Q. 
There are many possible deviations from the perfect 
model that could cause ZnCr 2 04 to order nonethe- 
less. These include further neighbor interactions and 
spin space anisotropy (|J|^|. Because the transition in 
ZnC^Cu is of the first order and involves a lattice dis- 
tortion, we suggest that finite lattice rigidity is an im- 
portant factor at the phase transition in this material. 
Consider the effect of tetragonal strain on magnetism in 
ZnCr204. It is well known that the exchange interaction 
between Cr 3+ ions whose oxygen coordination octahe- 
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dra share an edge is strongly dependent on the Cr-Cr 
spacing, r ]l9[| . Analysis of a series of chromium oxides 
indicates that dJ/dr ss 40 meV/A p0| . This implies that 
tetragonal strain e a > and e c < yields weaker AFM 
interactions in the basal plane, AJ± = r$e a dJ jdr = 0.06 
meV, and stronger AFM interactions between all other 
spin pairs A J|| = r (e Q + e c )/2dJ/dr — —0.04 meV. This 
asymmetry reduces the mean field energy of the ordered 
phase in ZnCr 2 4 by A(H S ) = (5AJ\\ - AJ x )/2 « -0.07 
meV/Cr relative to the mean field ground state energy 
in the cubic phase ^IJ. The result should be Neel order 
in tetragonal ZnCr204 below an ordering temperature 
that we denote Tpj t . Because order appears abruptly and 
simultaneously with the tetragonal strain we infer that 
T Nt > T c . This may be possible despite the modest value 
of A(H. S } because of the strong local constraints present 
when |T/0cw| << 1- In addition there could be other 
hitherto undetected lattice modifications at T c that also 
favor Neel order (see below). 

The magnitude of the lattice distortion is controlled by 
the need to balance the increase in lattice energy and the 
decrease in entropy against the decrease in the energy of 
the spin system. Equating the free energy F = (Tli + 
7i s ) — TS of the competing phases at T c implies that 



A(H S ) + A(Hi) - T C AS = 



(1) 



when cooling through T c . We can derive A(TC S ) from Fig. 
1 (b) and (c) using the first moment sum rule: [ ^2| 



A ^ = ~9^M • n in 

2 1 — sin Qr / Qr 



)AS{Q,uj)duj 
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Limiting the integral to 0.2 meV to 12 meV and aver- 
aging data for 1.3 A < Q < 2 A yields a value of 
A(H. S ) = —0.40(7) meV/Cr. From specific heat mea- 
surements H we find that AS = —0.107 i?ln4/Cr corre- 
sponding to — T C AS = 0.16 meV/Cr. The difference be- 
tween these numbers yields an estimate for the increase in 
lattice energy A(Hi) = 0.24(7) meV/Cr. The energy as- 
sociated with simple tetragonal strain J23| only accounts 
forv (c 11 (e 2 c +2e 2 a )/2+c 12 (e 2 a +2e a e c )) = 0.026 meV/Cr so 
there are likely additional modifications to the structure 
of ZnCr 2 04 below T c that help to stabilize Neel order. 

There are interesting analogies between the phase tran- 
sition in ZnCr204 and the spin-Peierls transition. In both 
cases the high T phase is near quantum critical and can 
lower its energy through a lattice distortion. In both 
cases the transition occurs from a strongly correlated 
paramagnet: T c <C Qcw- And in both cases low energy 
spectral weight is moved into a finite energy peak. There 
are also important differences that render the transition 
in ZnCr 2 04 a distinct new phenomenon in magnetism. 
The lattice distortion in ZnCr 2 04 drives the spin system 
into an ordered phase not a quantum disordered phase. 
The transition in ZnCr 2 C>4 is of the first order while the 
SP transition is of the second order. And the change in 



entropy at T c plays an important role in ZnCr 2 C>4, not 
in a SP transition. The central idea that finite lattice 
rigidity can drive a spin system away from quantum crit- 
icality however does carry over, and might be relevant for 
other frustrated magnets when symmetry breaking terms 
in the spin Hamiltonian fail to induce magnetic order. 
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